Numerical solution of partial differential equation